The dynamo effect - a dynamic renormalisation group approach 
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The Dynamo effect is used to describe the generation of magnetic fields in astrophysical objects. 
However, no rigorous derivation of the dynamo equation is available. We justify the form of the 
equation using an Operator Product Expansion (OPE) of the relevant fields. We also calculate the 
coefficients of the OPE series using a dynamic renormalisation group approach and discuss the time 
evolution of the initial conditions on the initial seed magnetic field. 
PACS no.: 47.25Cg,47.65.+a,91.25.Cw 

The existence of magnetic fields in the astrophysical objects has been an active area of research for quite sometime. 
This effect is known as the dynamo effect. Many mechanisms have been suggested to explain this. The dynamo effect 
' could be caused by a background turbulent fluid (see, for example [1]) or by a non-turbulent backgorund fluid motion 
£^ . (an interesting example of this type has been discussed in [2]). In this letter we, however, confine ourselves to the 
study of the first kind. 

In turbulent dynamo, to begin with one assumes that initial seed field is very weak i.e. most of the energy is 
contained in the kinetic part of the magnetised fluid. This immediately tells us that even though the velocity field 
. is influencing the time evolution of the magnetic field, the back reaction of the magnetic field on the velocity field 
is negligible as the ratio of the the Lorentz force to the inertial force is ~ B 2 /u 2 (this is an order of magnitude 
on : estimate of the two nonlinear term in the Navier-Stokes equation) where B and u are the magnetic and velocity fields. 

Basically, the velocity field being unaffected by the magnetic field, will time evolve according to the Navier-Stokes 
(NS) equation and magnetic field will time evolve according to the Induction equation, being influnced by the velocity 
field. However, this assumption of weak magnetic field will hold only during the initial transient. 

Conventionally, while deriving the dynamo equation one takes a two scale approach [1]: It is assumed that there 
are two scales of variations -(i) A global scale of variation L of 'mean' quantities and (ii)a scale of small wavelength 
fluctuations l . Mean quantities are defined to be quantities which are averaged over an intermediate scale a:l Q << 
a << L. In other words, fields have a small wavenumber ~ L _1 and a large wavenumber ~ l^ 1 components and the 
mean quantites can be calculated by averaging over the large wavenumbers. We split the velocity and the nagnetic 
£^ ' fields into mean and fluctuating parts: 

o ' 

O ■ U(x,f) = U Q (x,t) +u(x,t), < u >= (1) 

B(x,i) = B o (x,f)+b(x,i),<b>=0 (2) 

<> indicates an averaging over a, the intermediate scale. 
■ The time evolution of the magnetic fields is governed by the Induction equation [1,3]: 

- = Vx(uxB)|/iV 2 B (3) 

where y, is the magnetic viscosity. With the above definitons of the mean and fluctuating quantities, one can easily 
write down the equations for B D and b. The one for B Q is given by 



a 

T3 



dt 



V x (U x B ) + V x E + z/V 2 B Q (4) 



where E =< u x b >. It is easy to see that, even though < u > and < b > are zero individually, their product may 
yield a non-zero value because of a possible long range correlation between the two. Then a gradient expansion is 
performed over this electromotive force E. This has been justified by the argument that E and B Q are linearly related: 

dB Q1 d 2 B 01 

Ek = '^1-.; • i,— I -^-J-J^ + - (5) 

In other words, averaging u x b over the intermediate scale gives rise to a gradient expansion of the long wavelength 
part of the magnetic field. The first term is responsible for 'a-effect' and the second term for '/3-effect'. It can be 
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easily shown that when the background fluid motion is homogenous and isotropic, a is a pseudo scalar and f3 is a 
true scalar. The j3 term contributes to the effective viscosity ('turbulent diffusion'). Hence, including these a and (3 
effects, the equation for the large scale component of the magnetic fields becomes 

-^ = Vx(U o xB )|aVx B + t?V 2 B (6) 

where r\ — v + (3, the effective viscosity. This is the dynamo equation. 

Even though this equation has been very successful, the justification for the existence of a gradient expansion (which 
is a crucial part in the derivation) is not founded upon a strong basis. Here we attempt to give a rigorous basis to 
it: From a field theoretic point of view < u(ri) x b(ra) > (here u(ri) and b(r2) are the 'fluctuating' part of the 
total velocity and the magnetic field as defined in equations (1) and (2)) diverges as n — > r 2 . Hence, one can write 
down an Operator Product Expansion (OPE) [4] for the product of the fields Lt ri ^r 2 u ( r i) x b(r 2 ) an d consequently 
calculate Lt ri ^ r . 2 < u x b >. Since the LHS is linear in b, the RHS will have only odd powers of B (so that B — ► B 
property is retained). For the same reason, there cannot be any term having no B . Also, no term propertional to b 
can appear as b.u x b = (the LHS). However, no such condition is there for the velocity field as if u is a solution 
of the NS equation then — u is not. We write down the OPE as 

OB 

(u x b)i =Ei = a tJ B 01 + Pijk-g^- + - + [] (7) 

where [] indicates various composite operators permitted by the symmetry of the LHS. Since, expectation value of 
any composite operator is defined to be zero we have 

dB 

< u x b >= aijBoj + Piik-Q^- + •••• (8) 

which is in agreement with the form of the equation. 

In fc-space, Lt ri ^ r2 < u x b >= J q u(q) x b(k — q). Here, we take k to be small i.e. long wavelength limit and 
A > q > Ae _r with r positive, i.e. q belongs to some short wavelength band. The RHS becomes, J aij(q)B j(k — 
q) + J (3ijk(q)kjB ok (k — q) + ... (in principle, a, (3 etc. could be functions of fc). For calculating the OPE coefficients, 
we calculate the LHS by momentum shell integration, do a loop expansion and equate various powers of fc on both 
sides. We work in the incompressible limit. 

The Navier-Stokes equation is given by 

^ + (u.V)u = z,V 2 u-^P+f (9) 

at p 

V.u = (10) 

The last one is the incompressibility assumption. Here, v is the fluid viscosity, p the pressure and p the density. This 
in fc-space reduces to 

Ou a (k, t) + i-P a ^{k) J Uf 3l (q)u k (k -q) = -vk 2 u a + f a (11) 

where P Ql a 7 = Papk^ + P Q7 fc/3 and P Q( a is the projection operator. Since we are looking at turbulent dynamo, we take 
< f a (k,t)fp(k',t') >= L>i(fc)P Q(3 (5(fc + k')5(t - t') + D 2 {k)e af)l k 1 8{k + k')6(t - t'). D 1 and D 2 are even in fc. An 
explicit factor of i infront of D 2 term indicates that it is an odd parity breaking term and according to our previous 
analysis a-effect will be absent if D 2 = 0. 
The induction equation is given by 

This, in fc-space becomes 

db a (k, t) 
dt 



— = V x (u x b) + M V 2 b (12) 

= itaPjkpe^^x / u AJ (q)6 A (fc - q) - pk 2 b a (13) 

Jq 

At the tree level (fig.l), the diagram which will contribute is the following: We have, < (u x B) Q >=< 
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if qn e a01 e lS xe\r,T)(k-q) s [D 1 (q)Pi 3ri (q)+D2{q)ep vp q p ]Gl(k-q,u-n)CZ(q,uj)(k- where G b a (q, ft) = 

and C"(g, fi) = n i+ v i q i ■ We extract the a-term and the (3 term from the above expression in the following way: 
I. The a term: In the above integral, there is a term which is independent of k. This gives rise to the a term 

aB a (k) = J d\dne Mx e Xw e tinp {k ^ g) ^ _ fl) | ^ g)2 B T {k) (14) 

Now, ea/sje^sx^XrjT^prip — S a sS Tp + S ap 5 T s. Hence, one obtains 

47r 2 1 r 

a=——— D 2 (q) (15) 
3 v(v + ti) J q 

Since, the integral is infra-red divergent, we evaluate it by shell integration. 
77. The [3 term: The term proportional to k gives rise to the (3 term 

if3e\s T = i J d 3 qdnDi(q)e af3l e 7 s\^x V TPpr 1 (q)(k - q)s (16) 

Now, e Q/ 3 7 e 7 5A = 5 aS 5p\ - 5 a \8p\. Hence, 



4tt 2 1 f D 1 (q) 
3 v(v + ix) 



if 



Notice that a is pseudo-scalar and (3 is a true scalar. We see that even for thermal noise i.e D\ and D 2 being 
constants give nonvanishing a and (3. However, if we assume that the background fluid is fully developed turbulent, 
then Di(k) ~ fc~ 3 and D 2 (k) ~ fc -4 (this particular choice gives E(k) ~ fc~ 5 / 3 [5]). 

Having calculated a and /? from the tree level diagram, we now analyse the higher-loop diagrams. There are several 
distinct classes of them: (i) This type can be generated by decorating the velocity correlator or the response function. 
They generate corrections to the coefficients already obtained without decorations (fig 2a). (ii) This type of diagrams 
are generated by decorating the magnetic field response function. They give rise to non-zero values for the higher 
order coefficients of the OPE series as well as corrections to the coefficients already obtained without decorations (fig 
2b). These two classes of diagrams can be obtained by considering a 'self-consistent' tree level diagram. 

(iii) These type of diagrams appear by inserting two ub vertex (the advective vertex of the Induction equation) in 
the tree level diagram (fig 3a) or by inserting one uu and one ub vertex (fig 3b and 3c). A 1-loop diagram as shown 
in fig 3a is given as J G b a (k — q)C%(q)Ai(p) where 

A 1 {p) =G b (k-q- p)G b (k - p)C:(p) (18) 
Similarly, fig 3b and fig 3c are given by J q p G b (k - q)C^(q)A 2 (p) and J q p G b (k - q)C^(q)A 3 (q) where 

A 2 (p) = G b a (k - p)G u (p - q)C:{p) (19) 
MP) - G b a (k - p)G u {-p)C u {p - q) (20) 

We see that the sum of Ai,A 2 and A 3 are nothing but the 1-loop correction of the advective vertex of the Induction 
equation. Sum of A\,A 2 and A 3 are however not zero as we are calculating the vertex correction at finite external 
momenta. These diagrams give rise to 1-loop corrections to a and (3 (which are already present in the tree level) as 
well as the lowest order 7. It is easy to see that an n-loop diagram will give rise to terms upto the n-th term in the 
OPE series starting from a. 

A pertinent question at this stage is:How does initial correlations among the seed magnetic field (a dynamo is 
essentially an initial value problem) grow? Also, if the initial magnetic field correlations are parity symmetric, will 
it generate a parity breaking correlation after a finite time (possibly due to the parity breaking correlation of the 
forcing in the Navier-Stokes equation)? Consider a 1-loop diagram which would renormalise the initial magnetic field 
correlations: It is clear that inside such a 1-loop integral, one noise correlation arises from the NS equation and the 
second one is from the initial magnetic field correlation. Hence, if the N-S noise has a nonvanishing parity breaking 
term, it will give a non-zero contributuion to the 1-loop integral only if the initial seed magnetic field correlation 
also has a parity breaking term, as one requires two or even number 'epsilon' terms for the integral not to vanish. 
Consequently, we can say that unless the initial correlation has a parity breaking term, it will not be generated in the 
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future time. It is also clear that even if the intial magnetic field correlation is non-singular at t = 0, it will pick up a 
singular contribution at a later time (since the noise in the NS equation has a singular correlation). 

It is clear that due to the dynamo action, the initial magnetic field will grow. Now, once J B 2 / J v 2 becomes ~ 1, 
it is nolongcr justified to neglect the Lorentz force term, which we have neglegted so far. Adding this term to the N-S 
equation we have the following set of equations: 

5j + (u.V)u = ^V 2 u + +-^(VxB)xB + --+f (21) 
at vkp p 

V.u = (22) 
— = aV x B + V x (u x B) + tjV 2 B (23) 

V.B = (24) 



We choose 

< fi(k,t)fj(k!,1? >= [D^k- 3 + iD 2 e ijp k p ]S(k + k')8{t - t') (25) 



< B i (k,t = Q)B j (k' \t = 0) >= DP l:j k- 3 S(k + k') (26) 

It is clear the dynamo equation has two linear terms - one proportional to k and the other proportional to k 2 . 
Conventionally, one may like to keep all the linear term in the bare response function. However, that would give 
rise to a singularity at a finite real k. To avoid that problem, we keep only the viscosity term in the bare response 
function and treat the a term perturbativcly. Following Forster et al [6] we obtain the following recursion relations 
for the parameters: 

da 1 

— = [z-l + A 6 29 
dr olv(v + rjjA 1 -' 6 

One can see easily that a ~ fc -1 / 3 and r\ ~ fc~ 4 / 3 . Ai,...,As are numerical constants which determine the non- 
universal amplitudes. It is clear that all the exponents dynamic (=2/3) and roughness exponents are same as usual 
MHD turbulence. Only the nonuniversal amplitudes are different. Here rj is the effective viscosity (or the 'turbulent 
diffusion' as opposed to the bare molecular viscosity. This is the origin of the enhanced diffusion process in many 
astrophysical situations e.g accretion disk [7]. 

So, in conclusion, we have shown how the dynamo equation can be formally constructed using OPE. We have shown 
how to calculate the coefficients of the OPE series pertubatively in a loop expansion. In particular we have evaluated 
a and (3 in a momcntun shell cllimination method. We have also discussed renormalisation of the dynamo equation. 
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Figure Captions: Fig 1: Tree level diagrams for < u(q) x b(k — q) >. A solid line indicates a bare magnetic field 
response function, a broken line indicates a bare velocity response funtion, a 'o' joined by two broken lines indicates 
a bare velocity correlation function, a wavy line indicates a magnetic field, a solid triangle indicates a ub vertex. 

Fig 2a and 2b: 1-loop corrections to the previous tree level diagram which arise due to dressing of the bare velocity 
response/correlation functions and bare magnetic field response function, a 'X' indicates a uu vertex. Other symbols 
have same meaning. 

Fig 3a, 3b and 3c: 1-loop correction to tree level diagram arising out of 1-loop corrections to the ub vertex. Symbols 
have same meaning. 
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